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Quenched staggered light hadron spectroscopy from 48 3 x 64 at (3 = 6.5* 
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We report our light hadron mass calculation based on an increased statistics of 250 quenched gauge configu- 
rations on a 48^ x 64 lattice at j3 — 6.5. Quark propagators are calculated for each of these configurations with 
staggered wall source and point sink at quark mass values of m q — O.Of , 0.005, 0.0025 and O.OOf 25. We also did 
additional calculations to improve our understanding of systematic biases arising from autocorrelation, source 
size, and propagator calculations. Our earlier conclusions that the flavor symmetry breaking is reduced and the 
ratio rriN /m p (~ 1.25(4)) is small remains robust. 
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Understanding low energy properties of the 
strong interaction from first principles of quan- 
tum field theory is one of the main goals for 
the lattice quantum chromodynamics (QCD). 
Hadron spectrum is a typical example of such low- 
energy phenomena and precision determination of 
hadron spectrum can serve as a validity check of 
lattice QCD 0. 

Over the past few years, we have been calculat- 
ing light hadron masses using quenched approx- 
imation to lattice QCD on a large lattice vol- 
ume (48 3 x 64) with small lattice spacing (cou- 
pling constant of f3 = 6.5) and small quark mass 
(m, = 0.01, 0.005, 0.0025 and 0.00125) §. There 
have been numerous efforts (see for example |j],|| ) 
for quenched light hadron spectrum calculation, 
most of which required extrapolations with re- 
spect to lattice volume, lattice spacing, and quark 
mass to approach large physical hadron size, con- 
tinuum limit and light up and down quark mass. 
Controlling uncertainties arising from such ex- 
trapolations is difficult. In particular, extrapo- 
lations with regard to quark mass can be trou- 
blesome due to the fact that the chiral behavior 
of quenched theory is different from that of full 
theory 0. 

We try to reduce various systematic errors as- 



sociated with the extrapolations by calculating on 
a large lattice volume, small lattice spacing and 
small quark mass. Previous studies suggest that 
for quenched staggered spectrum, (3 — 6.5 lies in 
asymptotic scaling region. Thus the question lies 
on whether we can simulate large lattice volume 
at (3 = 6.5 and can reduce quark mass. 

We use a combination of Metropolis and over- 
relaxation sweeps for generating quenched gauge 
field configuration. As we noted last year, we 
increased the sweep separation between measure- 
ment from 1000 to 2000 sweeps. This was neces- 
sary for the lightest quark mass value of m q a = 
0.00125 to be free of autocorrelation, but not 
for the three heavier mass values. For hadron 
spectrum calculation, staggered quark wall source 



with 



0.01,0.005,0.0025 and 0.00125 and 
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point sink is used. For Dirac matrix inver- 
sion, standard conjugate gradient (CG) method 
is used. This year we tested the numerical robust- 
ness of hadron propagators by comparing output 
from a stricter (machine-accuracy) convergence 
condition. We found no significant difference. 

Table summarizes our results. The quoted 
errors are from \ 2 fit- The corresponding jack- 
knife errors are in agreement within 10 to 20 %. 
We measured nucleon masses from "corner wall" 
source and from "even point wall" source. How- 
ever since "corner wall" nucleon effective mass 
does not show noticeable plateau, we present nu- 



cleon mass from "even point wall" source. The 
Goldstone pion mass (it) and non-Goldstone pion 
mass (7T2 ) agree within their errors for all the four 
quark mass values. So do the two different rho 
mesons, p and p2 obtained as parity partners of 61 
and a\ respectively. The flavor symmetry break- 
ing in these channels is still not detectable despite 
the refined errors. In other words, our conclu- 
sion last year that the flavor symmetry breaking 
is not seen at (3 = 6.5 is robust and has been 
refined. The smaller errors are obtained partly 
because of the increased statistics: the wiggling 
of pion effective mass also reported last year has 
been reduced. In addition we have better con- 
trol in selecting plateaus in the effective mass 
plot. This became possible through our source 
size study (see Figure |l]). We clearly observe 




the effective mass in earlier time decreases as the 
source size is increased. Wall size dependence 
for m q — 0.01 pion effective mass is qualitatively 
the same. This behavior gives us a clear indi- 
cation of how to define a plateau and results in 
smaller and more reliable error estimate. Since 
the combination m^L takes the values of about 
2.9, 3.8, 5.3, and 7.6 respectively for quark mass 
m q a = 0.01, 0.005, 0.0025 and 0.00125, we do not 
have to worry much about finite- volume effect on 
the current lattice with L = 48 either. 

Figure [2] shows our Edinburgh plot. Without 
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Figure 2. Edinburgh plot at [3 — 6.5 for quark 
mass m q = 0.01, 0.005, 0.0025 and 0.00125. Nu- 
cleon masses from even-point wall source are used 
for m^/m p at each quark mass values. 



Figure 1. Nambu-Goldstone pion effective mass 
at /? = 6.5 on 48 3 x 64 lattice for m q = 0.00125. 
Three new different sizes for corner- wall, 12 3 (□), 
24 3 (+), 36 3 (O) are used in addition to the 

48 3 (x). 



that the effective mass for all the four source sizes 
eventually approach a common plateau. The un- 
wanted contribution from the excited states to 



any extrapolation with respect to quark mass, 
win / wip approaches experimental value when the 
lattice volume is large enough to admit light pion. 
Let us look for quenched chiral logarithm in 
our light hadron spectrum. In Figure |3j, we plot 
m^/rug as a function of m V2 . Instead of stay- 
ing flat, the data points seem to be rising, simi- 
larly to what has been observed previously. This 
could be a finite volume effect as suggested in 
ref. ||. Therefore, in Figure [|, we plot m% as 
a function of m„. The two curves are results 



Table 1 

hadron masses for m q a = 0.01,0.005,0.0025 and 0.00125 



particle 



0.01 



0.005 



0.0025 



0.00125 



7T 


0.1576(3) 


0.1114(3) 


0.0811(6) 


0.0610(7) 


n-2 


0.1570(5) 


0.1125(6) 


0.0836(10) 


0.0637(20) 


a 


0.321(3) 


0.320(6) 


0.340(10) 


0.291(10) 


P 


0.2430(7) 


0.222(1) 


0.216(2) 


0.209(4) 


P2 


0.2412(7) 


0.221(1) 


0.216(2) 


0.212(3) 


oi 


0.346(3) 


0.326(4) 


0.330(4) 


0.311(5) 


h 


0.344(4) 


0.326(6) 


0.333(8) 


0.349(16) 


N 


0.337(1) 


0.293(2) 


0.269(3) 


0.261(6) 


A 


0.397(3) 


0.382(3) 


0.368(5) 


0.360(7) 
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Figure 3. m^/m q (vertical axis, logarithmic scale) 
at f3 = 6.5 as a function of m- K2 (horizontal axis). 



from fitting to m 2 = cq + c\vn q following sugges- 
tions by H. The broken line is fitted result us- 
ing m q = 0.01, 0.005 and 0.0025 with x 2 /d.o.f ~ 
6.7. The dotted line is fitted result using m q = 
0.005,0.0025 and 0.00125 x 2 /d.o.f ~ 15.4. Both 
lines show non-zero intercepts (0.30(11) x 10~ 2 
for m q = 0.01,0.005,0.0025 and 0.80(10) x 
10~ 2 for m q = 0.005,0.0025,0.00125. However, 
the slopes change noticeably from 2.45(2) for 
m q = 0.01,0.005,0.0025 to 2.32(3) for m q = 
0.005, 0.0025, 0.00125. Since our pion mass is al- 
ready small for m q — 0.01, we think that the 
change in slope is less likely due to the neglected 
higher order terms, 0(m™)(n > 2). Thus, our 
data does not appear to agree with a finite vol- 
ume cutoff plus linear term picture of the pion 
mass. 



Figure 4. m 2 (vertical axis) at j3 = 6.5 as a func- 
tion of quark mass (horizontal axis) The lines are 
results from fitting to a form, m 2 = cq + c\m q 
with three data points only. 
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